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Abstract 

The effect of the recently obtained 2nd post-Newtonian corrections on the accuracy of estimation 

of parameters of the gravitational-wave signal from a coalescing binary is investigated. It is shown 
that addition of this correction degrades considerably the accuracj' of determination of individual 
masses of the members of the binary. However the chirp mass and the time parameter in the signal 
is still determined to a very good accuracy. The possibility of estimation of effects of other theories 
of gravity is investigated. The performance of the Newtonian filter is investigated and it is compared 
with performance of post-Newtonian search templates introduced recently. It is shown that both 
search templates can extract accurately useful information about the binary. 
PACS numbers: 04.30.+x,04.80.+x,95.85.5z,97.80.Af 
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1 Introduction 



It is currently believed that the gravitational waves that come from the final stages of the evolution of 
compact binaries just before their coalescence are very likely signals to be detected by long-arm laser 
interferometers 0|. The reason is that in the case of binary systems we can predict the gravitational 
waveform very well; and the amplitudes are reasonably high for sources at distances out to 200 Mpc. 
An estimate based on the number of compact binaries known in our galaxy and extrapolated to the rest 
of the Universe shows that there should be one neutron star compact binary coalescence per year out 
to the distance of 200Mpc |^. This estimate is a safe lower bound on the rate of binary coalescence. 
Arguments based on progenitor evolution scenarios suggest that there should be 100 of two neutron 
star coalescences, 5 neutron star - black hole coalescences, and 0.5 two black hole coalescences out to 
200Mpc Q. The waveform derived using the quadrupole formula has been known for quite some time 
A standard optimal method to detect the signal from a coalescing binary in a noisy data set and 
to estimate its parameters is to correlate the data with the filter matched to the signal and vary the 
parameters of the filter until the correlation is maximal. The parameters of the filter that maximize the 
correlation are estimators for the parameters of the signal. The detailed algorithms and the performance 
of the matched-filtering method in application to coalescing binary gravitational-wave signal has been 
investigated by several authors, e.g. [§[ 0, ||, |^, 0. It has recently been realized that the correlation 
is very sensitive even to very small variations of the phase of the filter because of the large number 
of cycles in the signal. Consequently the addition of small corrections to the phase of the signal due 
to the post-Newtonian effects decreases the correlation considerably. Thus the post-Newtonian effects 
in the coalescing binary waveform can be detected and estimated to a much higher accuracy than it 
was thought before This opens up new prospects but also considerable data analysis challenges 

for the LIGO, VIRGO, and GEO600 projects which are rapidly progressing. It was also found |^ 
that the post-Newtonian series is not converging rapidly for a binary near coalescence. Hence higher 
post-Newtonian corrections will affect the correlation. Currently three post-Newtonian corrections to the 
quadrupole formula are already known fl^ and the calculation of further ones is in progress. In this 
article we analyse the estimation of parameters of the 2nd post-Newtonian signal. This part of work 
complements a recent detailed analysis of the 3/2 post-Newtonian signal performed recently in Rcf . . 
We also examine the detectability of the post-Newtonian signal and estimation of its parameters using the 
Newtonian waveform as a filter. This filter can be used as the simplest search template. We compare the 
Newtonian search templates with the post-Newtonian search templates recently investigated in Ref . 10 . 

The paper is organized as follows. In the first part of Section 2 we present the gravitational wave signal 
from a binary system to the currently known 2nd post-Newtonian order. In this work we analyse the signal 
in the "restricted" post-Newtonian approximation (i.e. only the phase of the signal is given to the 2nd 
post-Newtonian accuracy whereas the amplitude of the signal is calculated from the quadrupole formula) , 
we assume circularized orbits, and we assume the spin parameters to be constant. In the second part we 
briefly describe the optimal method of detection of such a signal in noise and the maximum likelihood 
(ML) method to estimate the parameters of the signal. We derive a number of properties of the ML 
estimators of the parameters of our signal and we examine the bounds on their variances. Our analysis is 
based on the Cramer-Rao bound. In the third part we give the approximate rms errors of the estimators 
for the signal at various post-Newtonian orders. In the fourth part of Section 2 we consider the effects 
of other theories of gravity and their detectability from gravitational-wave measurements. We consider 
Jordan-Fiertz-Brans-Dicke theory and Damour-Esposito-Farese biscalar tensor theory. In Section 3 we 
consider the so called "search templates" introduced in Ref.[Q. These are simple filters containing as 
few parameters as possible to effectively detect the multi-parameter signals. In the first part of Section 3 
we analyse the simplest search template - the Newtonian filter which is the waveform of the gravitational 
signal from a binary in the quadrupole approximation. We examine the Newtonian filter as a tool both 
to detect the signal and also to determine its nature. In the second part of Section 3 we compare the 
Newtonian filter with the other search template analysed recently based on the full post-Newtonian 
signal. In Section 4 we summarize conclusions from our results. A number of results is left to appendices. 
In Appendix A we examine the first order effects on the phase of the signal due to eccentricity. In 
Appendix B we give numerical values of the covariance matrices at various post-Newtonian orders. In 
Appendix C we give certain detailed formulae for the Damour-Esposito-Farese theory. In Appendix D 
we briefly review the theory of optimal detection of known signal in noise and we generalize it to non- 
optimal detection. In Appendix E we give a useful analytic approximation to the correlation integral of 
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the optimal filter with the signal from a binary. 
The units are chosen such that G = c = 1. 

2 Post-Newtonian effects 

2.1 Gravitational wave signal from a coalescing binary 

Let us first give the formula for the gravitational waveform of a binary with the three currently known 
post-Newtonian corrections. We make the following approximations. We work within the so called 
"restricted" post-Newtonian approximation i.e. we only include the post-Newtonian corrections to the 
phase of the signal keeping the amplitude in its Newtonian form; this is because the effect of the phase on 
the correlation is dominant. The inclusion of post-Newtonian effects in amplitudes will not qualitatively 
change our results. Due to the effect of rapid circularization of the orbit by radiation reaction one can 
assume that the orbit is quasicircular. For example in the case of the gravitational wave signal from the 
Hulse- Taylor binary pulsar at the characteristic frequency of the detector for such a signal of around 47Hz 
the eccentricity e would be ~ 10~^. Moreover the first order contribution to the phase of the signal due to 
eccentricity goes like e^. Nevertheless for completeness we include the first order correction to the signal 
due to eccentricity in our formulae. We give a detailed derivation of this correction in Appendix A. We 
neglect the tidal effects. All tidal contributions to the gravitational wave signal from a coalescing binary 
were estimated to be small [^2[ |5) . There is also a small additional contribution to the phase due to tail 
effects which detectability has been considered in detail and was found to be small. This correction 
is formally of the 4th post-Newtonian order and consequently we neglect it in the present analysis. 

With these approximations the waveform, as a function of time, is given by the following expression, 

hit) = A f(tf/^ cos[2^ j'^ f{t')dt' ^ </.], (1) 

where 



(2) 

5 R 

and where 4> is an arbitrary phase, /i and m are the reduced and the total mass of the binary, respectively; 
ta is a time parameter and R is the distance to the source. A is the rms average amplitude over all Euler 
angles determining the position of the binary on the sky and the inclination angle between the plane of the 
orbit of the binary and the line of sight. The rms amplitude A is 2/5 of the maximum possible amplitude. 
The characteristic time for the evolution of the binary to the currently known 2nd post-Newtonian order 
is given by 

_ ./ _ 5 1 1 
^''''^ ■ df/dt 96 ^m2/3 (^/)8/3 ^•'> 

[1 _ + {— + — -^^)(7rm/)2/3 - (47r - 5o)(W) + (4) 

I 24 / 19/9 ^336 4 TO ^ ojy jjT vy 



,3058673 5429 /1 617, /i,, to,, 

( 1 — H ( — ) + Ss — )Umfr'^] 

^016064 1008 TO 144 TO^ ^'^ ^ 

where is the asymptotic eccentricity invariant, 

eo is the eccentricity of the binary at gravitational frequency /o (see Appendix A for derivation and 
explanations). The quantities Sq and Sg are spin-orbit and spin-spin parameters respectively. They are 
given by the formula 

113 25 TO2 TOl 

So = -TTrisi + S2) + -r(si + S2 — ) (6) 

12 4 TOl TO2 

247 721 
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Si = — , S2 = — (8) 

si = L-Si, S2=L-S2, (9) 

where L is the total orbital angular momentum and Si, S2 are the spin angular momenta of the two 
bodies. The terms in square brackets in Eq.(^ are respectively: at lowest order, Newtonian (quadrupole); 
at order f~^^^^, lowest order contribution due to eccentricity (see Appendix A); at order f"^^^, IPN JTzt ; 
at order /, the non-linear effect of "tails" of the wave (47r term) |3[ |2l|), and spin-orbit effects 
H; and at order /"/^ 2PN H and spin-spin effects [||. 

In general the spin parameters vary with time. It was shown that So is nearly conserved, it never 
deviates from its average value by more than ~ 0.25. Moreover the time dependent part of the spin 
parameter is oscillatory what reduces considerably its influence on the phase of the signal |l^. In this 
work we shall assume that both spin-orbit and spin-spin parameters are constant. We also neglect the 
effect of the precession of the orbital plane due to spin on the waveform. The effects of the spin on the 
waveform of the signal from an inspiralling binary have been investigated in detail in ||2^ . If we take the 
available estimate of the moment of inertia for the pulsar in the Hulse- Taylor binary and assume masses 
of the neutron stars in the binary of 1.4 solar masses then So — 4.8 x 10^^ and Sg — 2.4 x 10~^. If such 
values are typical then spin effects will make negligible contributions to the phase of the signal. However 
this may not be the case for binaries containing black holes. Moreover if cosmic censorship is violated and 
black holes rotate at a higher rate than allowed by maximally rotating Kerr black hole the spin effects 
will significantly affect the gravitational waveform. 

In the analysis of the detection of the above signal and estimation of its parameters it is convenient to 
work in the Fourier domain. The expression for the Fourier transform of our signal in the stationary 
phase approximation is given by (cf.[|l], ^, [lOt ) 

/i = i/"^/^expi[27r/ta - 0- 7r/4+ (10) 
5 

^(a(/;/a)fc + ae(/;/a)fce +ai(/;/a)fcl + a-i/2{f\fa)h/2+a,2{f;fa)k2)], 

for / > and by the complex conjugate of the above expression for / < where 

1 1 /ii/2mi/3 

^ = /onM 19. T, ' (11) 



(30)1/2 ^2/3 ^ 

273 ' 



^--ir, (12) 



= ^^e2(7^/o)l9/^ (13) 

1J43 11 , 
/i 336 4 

^3/2 = (47r-So), (15) 

to4/3 / 3058673 5429 617 . 
^^ = ^(,1016064 + 1008"+ 144" +7j' 
. ^ 9 1 3 ^/ 3 1 

40 (^/)5/3 + 8 (7r/,)s/3 5 (^/,)5/3 ' ^ 
... X 157 / 81 1 9 vr/ 9 1 \ 

«e(/; /a) - 1^1462 (^/)34/9 + 43 (7r/,)43/9 ~ 34 (^/a)34/9 j ' ^ 



ai(/;/a) = + ^7^^ ^, (19) 



^V2U-Ja)- (^0(^^)2/3+5(^/^5/3 2(^/,)V3)' (^0) 

9 1 3 tt/ 3 
"^(^'^''^ - aJ^^ + 4 (7^ (^^i^ ^^^^ 



hold and where fa = f{ta)- 
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The stationary phase approximation, Eq.(pl|), is an excellent approximation of the Fourier transform of 
the signal for frequencies which are not influenced by the finite time window of the measurement. In the 
above expressions for the gravitational wave signal from a binary we can make an arbitrary choice of the 
time parameter and the phase of the signal. 

We also point out that going from the time to the frequency domain we have made yet another approxi- 
mation. Namely we have taken the modulus | ft, | of the Fourier transform to be the Newtonian one i.e. 

I h |~ In the stationary phase approximation | h \ goes like j and consequently by Eq.(153) 

there would be other powers of frequency due to the post-Newtonian effects. We neglect those additional 
terms since the post-Newtonian corrections to the phase have the dominant effect. The inclusion of the 
post-Newtonian amplitudes to the signal will not qualitatively change the results of this work. 

A convenient parameter is the chirf mass defined as Ai := k~^^^. In the quadrupole approximation the 
gravitational wave signal from a binary is entirely determined by the chirp mass. 

We shall consider three models of binaries: neutron star/neutron star (NS-NS), neutron star/black hole 
(NS-BH), and black hole/black hole (BH-BH) binaries with parameters summarized in Table I. 

Table I. Numerical values of the parameters of the three fiducial binary systems. Black holes are of 10 
solar masses and neutron stars are of I.4 solar masses. Spin parameters are assumed to be constant. 
Spin for neutron stars was calculated from the typical estimate of the moment of inertia I for a neutron 
star of I ~ lO'^^ kg rr? . 



Binary 


mi[M0] 


m2[Mo] 


M[Mq] 


Sl 


S2 


So 




1. NS-NS 


1.4 


1.4 


1.2 


1.5x10-^ 


1.5x10"^ 


4.8x10-^ 


2.4x10-5 


2. NS-BH 


1.4 


10 


3.0 


3.0x10-5 


0.38 


4.0 


3.5 xlO-'' 


3. BH-BH 


10 


10 


8.7 


0.13 


0.13 


3.9 


0.15 



B 



Binary 


k[M-"''] 






k2WC] 


1. NS-NS 


0.72 


4.1 


25 


13 


2. NS-BH 


0.16 


2.0 


16 


15 


3. BH-BH 


2.7x10-2 


0.58 


4.7 


7.7 



where Mq means solar mass. 

For neutron stars we calculated the spin using the available estimate of the moment of inertia for the 
neutron star in the binary pulsar PSR1916-f 19. We have taken black holes to be spinning at half the 
maximum rate (i.e. Si = 0.5m^ /m^). The orbital momenta vectors were assumed to be parallel to the 
spin vectors. 

To have an idea of the size of the post-Newtonian corrections in the gravitational wave signal from a binary 
when it enters the observation window of the laser interferometer we have evaluated the characteristic 
time T2PN for the above three models at the frequency = 47 Hz which is the characteristic frequency 
of the detector for this signal (see below). We have made explicit the contributions to the characteristic 
time from the three post-Newtonian corrections. 

rjpjv = 44(1 + 0.046[from Ipn] - 0.025[from 3/2pn] + 0.0012[from 2pn])sec (22) 
t|p^ = 9.9(1 -I- 0.10[from Ipn] - 0.071[from 3/2pn] + 0.0060[from 2pn])sec (23) 
■^2^ = 1-7(1 + 0.17[from Ipn] - 0.12[from 3/2pn] + 0.018[from 2pn])sec (24) 

One concludes from the above numbers that for the earth-based laser interferometers post-Newtonian 
corrections are significant. Moreover several things are apparent. The quadrupole term is dominant for 
all the three models. This indicates a very good accuracy of the quadrupole formula even in the regime 
of strongly gravitating bodies. This has been noticed in other studies for example in the numerical 
investigation of the gravitational wave emission from the two black hole collisions H]. The difference 
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in size between the 1st post- Newtonian correction and the 3/2 post-Newtonian correction (tail term) 
is rather small. They differ by a factor of 2 for NS-NS binary and only by a factor of around 1.5 for 
binaries with a black hole. The second post-Newtonian correction is noticably smaller then the 3/2 post- 
Newtonian correction. The difference varies form a factor of 20 for a NS-NS binary to a factor of 7 for a 
BH-BH binary. The convergence of the post-Newtonian series appears to be worst for BH-BH binaries 
and in this case it would be desirable to have accurate numerical waveforms and not only the ones based 
on the post-Newtonian approximation. Such waveforms should be available as a result of the numerical 
projects such as Grand Challenge project currently under way in the United States. 



2.2 Detection of the signal and estimation of its parameters 

For the purpose of this investigation we shall use a fit to the total spectral density Sh {f) of the noise in 
the advanced LIGO detectors, devised in jl^. This fit comprises seismic, thermal, shot, and quantum 
noises in the detector. 

ShW = SoiUo/ff + 2(1 + (///o)'))/5, (25) 

where /o = 70Hz and 6*0 = 3 x lO^'^^Hz^^. It is an excellent approximation to the detailed formulae 
for various noises given in The sensitivity function Sen{f) of the detector is defined as \/Sh{f)- 
The sensitivity function has the maximum at frequency fo given above and its half width half magnitude 
(HWHM) (Jo is around 48Hz. 

To determine whether or not there is a signal in a noisy data set we use the Neyman-Pearson test (see 
Appendix D). When the noise in the detector is Gaussian the Neyman-Pearson test is the correlator test. 
It consists of linear filtering the data with the filter which Fourier transform is the Fourier transform 
of the signal divided by the spectral density of the noise The signal-to- noise ratio d that can be 

achieved by optimal filtering is given by d = (/i|ft,)^/^ where following jl^ the scalar product {hi\h2) is 
defined by 

{hl\h2)=ml ^^df, (26) 



where 3? denotes the real part. Thus we have 



We shall call the integrand of the above signal-to-noise integral signal sensitivity function and we denote 
it by Ind{f). This function has the maximum at the frequency where — 47Hz and its HWHM a'^ 
is ~ 26Hz, around half of that of the sensitivity function. This is the signal-to-noise ratio after filtering 
of the data. We see that linear filtering introduces an effective narrowing of the detector bandwidth |Q . 

In the case of our chirp signal the linear filtering increases the signal-to-noise ratio by an amount given 
roughly by the square root of the number n{f) of cycles spent near the frequency where n(f) is defined 

by i 

Consequently the effectiveness of matched filtering falls with the chirp mass. On the other hand the 
amplitude A of the signal increases with the chirp mass like Ai^^"^ and the overall factor in the signal-to- 
noise ratio increases as Ai^^^. This is born out by the amplitude A of the Fourier transform. Thus the 
probability of detection of binaries with the same rate of occurrence increases with the chirp mass. 

To estimate the parameters of the signal it is proposed to use the maximum likelihood estimation (MLE) 
p8| . It is by no means guaranteed that this is the best or the ultimate method. It may sometimes fail to 
give an estimate and other methods may lead to more accurate estimates. The MLE method consists of 
maximizing the likelihood ratio with respect to the parameters of the filter. In the case of the Gaussian 
noise the logarithm of the likelihood ratio A is given by ||2^ 

\nA^{x\hF)~^{hF\hF) (29) 
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where hp which we call the filter has the form of the signal but with arbitrary parameters and x are the 
data. We assume that the noise n in the detector is additive i.e. x = h + n. The maximum likelihood 
(ML) estimators of the parameters of the signal are given by the following set of differential equations 
providing that one can differentiate under the integration sign of the scalar product defined above. 

{x-hF\hF,i)^0, (30) 

where is the derivative of hp with respect to the ith parameter. Rarely these equations can be 

solved analytically. It was shown that in the case of the signal from a binary within the stationary 
phase approximation analytic expressions can be obtained for the maximum likelihood estimators of the 
amplitude and the phase. 

The ML estimators are random variables since they depend on the noise. It is important to know the 
statistical properties of these estimators and their probability distributions so that we can determine how 
well they estimate the true values of the parameters. The most important quantities are the expectation 
value of the estimator and its variance. We would like to have the expectation value of the estimator to be 
as close as possible to the true value of the parameter and we would like the variance of the estimator to 
be as small as possible. The difference between the expectation value of an estimator of a parameter and 
the true value of the parameter is called the 6ms of the estimator. The ML estimator is not guaranteed 
to be either unbiased or minimum variance. We have the following useful general inequality called the 
Cramer-Rao inequality ||2^ that gives lower bound of the variance of estimators. Let (6*^) be a set of n 
parameters and let 9j be one of the parameters then the variance of its estimator 9j satisfies the following 
inequality 

Var[ei] > (T-'^)ija^a^, (31) 

where a' and F*-' are given by 

where E is the expectation value. The matrix T is called the Fisher information matrix and its inverse 
is called the covariance matrix. One easily sees from the above inequality that when an estimator 9i is 
unbiased then the lower bound on its variance is given by the (II) component of the covariance matrix. 
For this inequality to hold certain mathematical assumption must be fulfilled p9[ |. 

1. The likelihood ratio must be a differentiable function with respect to all the parameters 9i. 

2. The order of differentiation with respect to parameters and the integration in the expectation value 
integral must be interchangeable. 

3. The variances of the estimators must be bounded. 

4. The Fisher information matrix must be positive definite. 



The Cramer-Rao inequality is very general. It holds no matter what is the probability distribution of the 
data and it applies to any estimator providing the regularity conditions mentioned above are fulfilled. 
The above inequality guarantees only that the variance of an estimator is greater then a certain amount. 
It is important for us to know how well the right hand side of the Cramer-Rao inequality approximates 
the actual variance of an estimator. It was shown ^ that in the case of Gaussian noise and in the 
limit of high signal to noise ratio d to the first order the maximum likelihood estimators are Gaussian and 
moreover they are unbiased and their covariances are given by the covariance matrix defined above. In 
statistical literature there also exists a series of refined Cramer-Rao bounds called Battacharyya bounds 
p9| . However in our case a useful approach to have an idea of the accuracy of the Cramer- Rao lower 
bound is given in Ref. [lO| ] where the maximum likelihood equations were solved iteratively and a formula 
for the covariance matrix of the ML estimators was derived to one higher order then given by the inverse 
of the Fisher information matrix. This formula can be treated as an approximation to the variances of 
the ML estimators by a series in 1/d where d is the signal-to-noise ratio. The first order terms given 
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by the inverse of Fisher matrix go as l/d? and the correction terms go Hke l/d'*. Consequently one can 
expect that for signal-to-noise ratios of 10 or so the diagonal elements of the inverse of the Fisher matrix 
give variances of the ML estimators to an accuracy of few %. 

We shall show that the set of parameters that we have chosen for our chirp signal has particularly useful 
properties. Note that the phase of the Fourier transform is linear in the phase, the time parameter, 
and the mass parameters ki. We shall call these parameters phase parameters. Moreover the Fourier 
transform is linear in the amplitude parameter A. The maximum likelihood estimators are those values 
of the parameters that maximize the likelihood ratio. The expectation value of the log likelihood is given 

by 

E[\nK\^{h\hp)-]^{hF\hF), (34) 

where is called the correlation function and is denoted by H. Using the stationary phase approx- 

imation to the Fourier transform of the signal H is given by the integral 

H{At, A0, Afc, Afce, Afci, AA:3/2, Afca) = (35) 

^^j^7/3 cos[2^/At - A(/. -(a(/; /a) Afc + a^if; fa)Ake 

ai(/; fa)Aki + a3/2(/; fa)Aky2 + a2(/; /a)Afc2)], 

where At means the difference in time parameters of the signal and the filter. The expectation of 
the log likelihood ratio depends on the phase parameters only through the correlation integral since 
{hp\hp) — i?(0, 0, 0, 0, 0, 0, 0) = (P where d is the signal-to-noise ratio. We see that the correlation 
function depends only on the differences between the values of the phase parameters in the signal and 
the filter and it has the maximum when the differences are zero. Moreover the value of the correlation is 
the same if we move by the same amount from the maximum in any direction for a given parameter i.e. 
H{—At, 0, 0, 0, 0, 0, 0) = H{At, 0, 0, 0, 0, 0, 0) and so on for all phase parameter^ This property means 
that the probability distribution of any estimator of the phase parameter will be an even function of the 
difference between the estimator and its true value. In other words the probability distributions of the 
estimators of the phase parameters are symmetric about their true values. Consequently we have 

mi E[{§i - OiY] = for 1 odd (36) 

and moreover for 1 even the moments mi are independent of the true values of the phase parameters. 
Thus the ML estimators of the phase parameters are unbiased (this is immediate from Eq.^ for / = 1) 
and the covariance matrix of the estimators of the phase parameters is independent of their values. The 
probability distributions of the phase parameters will depend on the signal-to-noisc ratio. We know that 
for large signal-to-noise ratio they will tend to Gaussian probability distributions. The estimator of the 
amplitude parameter is biased nevertheless by the symmetry property of the probability distributions of 
the phase parameters its bias is independent of the values of the phase parameters. These properties of the 
parameters can be also seen explicitly from the first two terms of the series solution of the ML equations 
(Eq.^o|) given in rcf. [0. The properties of our chosen set of parameters greatly simplify calculation of 
the Cramer-Rao bounds. In our case the Fisher information matrix F is given by 

f) J-f 

where S refers to the parameters of the signal and F to the parameters of the filter. The inverse of the T 
matrix is called the covariance matrix and is denoted by C. It is easily seen that F'^* components are all 
equal to zero when i A. Thus the amplitude parameter decouples from the phase parameters. Because 
the phase parameters are unbiased the lower bounds of their variances are given just by the appropriate 
diagonal elements of the covariance matrix C. In the case of the amplitude parameter the Cramer- Rao 
bound is given by VarA > b'{A)/r^^ where b'{A) is the derivative of the bias of amplitude parameter 
w.r.t. amplitude and F"^"^ = (P /A^. Note that F'^'^ is independent of A. This is a consequence of the 
linearity of the signal in the amplitude. 

It is clear from the linearity of the function H in the differences A6 that the F matrix is independent of the 
values of the phase parameters. Thus the Cramer-Rao bound on these parameters is also independent of 

^We are indebted to Dr. J. A. Lobo for this observation, see also ||2f3| p. 276. 
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the values of the parameters. From the argument above we know that this holds not only for the bounds 
on the variances but also for the variances themselves. 

To obtain the maximum of the correlation each phase parameter of the filter has to match a corresponding 
parameter in the signal (see Eq.|36|). Thus by linear filtering we shall get estimates of the time parameter 
ta, phase, and the mass parameters ki. In the filter one can always make an arbitrary choice of the time 
parameter ta- For example instead of choosing ta as the time at which frequency is fa one can choose 
time ijj as the time at which the frequency is equal to This new choice is equivalent to the following 
transformation 

t'a=ta+ Saka + Seke + Siki + (53/2^3/2 + ^2^2, (38) 
0' = + 6oko + 6eke + 6iki + (53/2^3/2 + (^2^2, (39) 

where 

- 5 1 1 

= 256^(^i^" K^^' ^^^^ 

^° ^ 16^(^i^"W^^' ^^^^ 
. ^ 785 1 

~ 11008 ^(7r/a)43/9 (7r/^)43/9^' ^ ^ 

4352^(7r/,)34/9 (^/^)34/9^' ^4^^ 

= Ife^(^-(^^' ^''^ 

Si - ^(^-r^n)' (45) 

48 Tlfa {Tlf'a) 

^ 32^(^i^ " (^i^^' ^^^^ 

- _ 5 1 1 

= m^w^" (^T^^' (*^^ 

5 1 1 
= l6((W^~W^^" ^^^^ 

The mass parameter frequency functions ai{f;fa),{i — 0,1,3/2,2) in Eq.^ are then transformed to 
Oj(/;/a)- The mass parameters remain invariant under the above transformations. By linear filtering 
with the template parametrized by the new time parameter and the new phase given by the above 
transformation we estimate the new time parameter t'^ and the new phase (j)' but the same mass parameters 

ki . 

There is also a particularly simple parametrization of the signal. Let us rewrite the Fourier transform of 
the gravitational wave signal from a binary in the following form 

h = Af~''/^expi[2nftc- (t)c~TT/4:+ (50) 
3 k 4239 ke ^ 5 ki 3 A:3/2 ^15 ^2 ^ 



128(7r/)5/3 11696 (7r/)34/9 QQnf 32 (7r/)2/3 64 (7r/)i/3J' 

(for / > and by the complex conjugate of the above expression for / < 0) where tc and (j)c are coalescence 
time and phase respectively and they are given by 

5 _fc 785 fee 5 fci 1 fc3/2 5 k2 

' " 256 {nfa f/^ 110008 (^/)43/9 + 192 {irfa)^ 32 {nfa)^/^ 128 (7r/a)4/3 ' ^ ' 

^ _ ^ , 1 ^ '^^^ ^"^ , 5 /Ci 5 fc3/2 5 fc2 



16(7r/a)5/3 4352 (7r/)34/9 48 7r/a 32 (7r/a)2/3 16(7r/,)i/3- 
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Coalescence time amd coalescence phase are obtained when the time parameter t'^ is such that the 
corresponding frequency is infinite which occurs when the two point masses coalesce. We can estimate 
the coalescence time and the coalescence phase of the template if we filter for combinations of the time 
and phase parameters with the mass parameter given precisely by right hand sides of Eqs.(^l|) and (p^. 
There is also a transformation of the phase that we shall find useful (see next Section). 

0" = - 2T:fmta. (53) 

where fm is some arbitrary constant frequency. Using the new phase parameter in the filter given by above 
transformation we shall estimate a new value of the phase shifted by the amount 27Tfmta- It is not difficult 
to show that all the above transformations do not change the CR bound on the mass parameters however 
the transformation Eq.(^) changes the bound for time and phase parameters whereas transformation 
Eq.(^3|) changes the bound on the phase. We can use the freedom of these transformations in the filter 
to obtain better accuracies of estimation of the time and the phase parameters. 

2.3 Numerical analysis of the rms errors of the estimators 

First of all we investigate the influence of the increasing number of post- Newtonian parameters on the 
accuracy of their estimation. To this end we have calculated the covariance matrices for the signal con- 
taining only the quadrupole term, then covariance matrices for 1st post-Newtonian, 3/2 post-Newtonian, 
and 2nd post-Newtonian signal, and finally for the 2nd post-Newtonian signal with first order contribution 
due to eccentricity. The results are summarized in Table II where we have given rms errors of the phase 
parameters. We have given the rms errors for the time and phase of coalescence tc and (j)c respectively. 
We have also determined the frequency fmin for which the error in the time parameter is minimum and 
we have given the minimum error Atmin in the time parameter and the corresponding error A(j) in phase. 
We have considered a reference binary of = lAf© located at the distance of lOOMpc. We have taken 
the range of integration from lOHz to infinity. The signal-to-noise ratio for such a binary is around 25. 

Table II. The rms errors for the phase parameters at various post-Newtonian orders for a reference 
binary of chirp mass of 1 solar mass at the distance of lOOMpc. Expected advanced LIGO noise spectral 
density is assumed and the integration range from lOHz to infinity is taken giving signal-to-noise ratio of 
around 25. 



Atm [msec] 




Ate [msec] 


A<t,c 


Ak[M~^^^] 


Afci[M-i] 


Ak3^^[M-''^''] 


Afc2[Af-i/^] 




0.14 


0.073 


0.17 


0.10 


8.3 xlO-" 










0.15 


0.087 


0.27 


0.33 


4.0 xlO-'^ 


5.8 xlO-3 








0.18 


0.14 


0.54 


1.9 


1.7 xlO"" 


0.70 xlO-i 


0.52 






0.24 


0.14 


1.6 


24 


6.6 xlO-i 


0.50 


7.2 


28 




0.25 


0.17 


2.3 


45 


2.3 xlO-3 


1.3 


17 


59 


1.2 xlO"*^ 



The above bounds scale exactly as the inverse of the signal-to-noise ratio d and they do not depend on 
the numerical values of the phase parameters. From the above table we see that increasing the number 
of post-Newtonian corrections and parameters we filter for decreases the accuracy of estimation of the 
parameters independently of the size of the post-Newtonian correction. Thus searching for a negligible 
correction due to eccentricity increases the rms error in other parameters by over 100%. 

For completeness in Appendix B we give the numerical values of covariance matrices for the phase 
parameters at various post-Newtonian orders and the corresponding values of the frequency fmin- 

As we have indicated above the estimator of the amplitude parameter is biased however if one takes the 
expansion of the variance of the estimator in the inverse powers of the signal-to-noise ratio (see |lQ] for 
a general formula) then the leading term for the variance of the amplitude is just l/T'^'^ where T^'^ is 
independent of A. The higher order corrections to the CR bounds of the amplitude go like l/d^ and they 
do depend on the value of the amplitude. As an amplitude parameter we find convenient to choose Aq 
given by 

= (54) 

flOOMpc 
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where Mq is the chirp mass m the units of solar masses and ^looMpc is the distance in the miits of 
lOOMpc. For our reference binary the ampUtude = 1 and thus the approximate rms error in its ML 
estimator is A^/d ~ 1/25 = 0.04 and as explained above this last number is independent of the true 
value of the amplitude. 

It is important to assess the accuracy of estimation of the physical parameters of the binary, i.e. the 
two masses of its members and the spin parameters So and Sg- This means that we have to make a 
transformation to a different parameter set. A nice property of the ML estimators is the following. Let 
9i be the maximum likelihood estimators of the set of parameters 9i. Let f{Oi) be a function of the 
parameters then f{9i) is the maximum likelihood estimator of the function / (see Ref[^). However 
it is not true in general that if estimators of the old parameters are unbiased then the new parameter 
is unbiased as well. Consequently by transforming the bounds of the old parameters one will not get 
the Cramer-Rao bound on the new set of parameters. However we know that Cramer-Rao bounds are 
approximately equal to the true variances in the limit of high signal-to-noise ratio d, correction terms 
being of the order of l/d^. Hence by transforming the C-R bounds one gets the rms errors of the 
estimators accurate to the order 1/d. Another important point is that the transformation to the new 
parameter set may be singular. Then the determinant of the F' matrix for the new set of parameters is 
zero and thus T' is not positive definite, consequently the Cramer-Rao inequality does not hold. A way 
to get errors of estimators of the new parameters in such a case could be to attempt to calculate the bias 
and the variance directly from some approximate probability distributions for the estimators (see ref . ||l^ 
for such treatment to determine the accuracy of the distance to the binary). It may happen however 
that the probability density function is such that the expectation value and the variance do not exist (an 
example is Cauchy probability distribution) and then one may have to use another measure of bias and 
error, e.g. median and interquartile distance. The other method proposed in [ p^ is to use confidence 
intervals. We shall return to this problem in the future work pT[ . 

The transformation from the 4 mass parameters fc/ to new parameters - total mass (m) , reduced mass (/x) 
and the spin parameters So and Sg is regular. Thus we can obtain approximate values of the errors of the 
estimators of the reduced mass, the total mass and the spin parameters. However the transformation from 
m and /U to individual masses mi and m2 is singular (determinant of the Jacobian of the transformation 
is zero when masses are equal, see Ref. [0 ) . Consequently the errors in the determination of the masses 
cannot be obtained from the C-R bounds calculated above. 

In Table HI we show the degradation of the accuracy of estimation of the chirp mass, the reduced mass, 
and the total mass with the increasing number of parameters in the signal for the NS-NS binary at a 
distance of 200Mpc. 

Table III. Degradation of the accuracy of estimation of the chirp mass, the reduced mass, and the total 
mass for the fiducial neutron star - neutron star binary with increasing number of parameters in the 
signal. 2 pN means that the phase of signal is taken to 2nd post- Newtonian order with spin parameters 
included and we maximize the correlation of the signal with a template matched to the signal for all the 
phase parameters. 



pN order 


AM/M 


Afi/fi 


Am/m 


1 pN 
3/2 pN 

2 pN 


0.0054% 
0.023% 
0.080% 


0.55% 
6.4% 
42% 


0.81% 
9.6% 
63% 



For the calculation of the numbers in the table above and all other tables in the remaining part of this 
Section we have taken the range of integration in the Fisher matrix integrals to be from lOHz to the 
frequency / = (6^/^7rm)~^ corresponding to the last stable orbit of the test particle in Schwartzschild 
space-time. This may very roughly correspond to the last stable orbit in a binary |^ , |33| . 

In Table IV we give the signal-to-noise ratios and the Cramer-Rao bounds for the mass and the spin 
parameters in percents of their true values for the 2nd post-Newtonian signal for our three representative 
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binary systems at the distance of 200Mpc. We have also given the improvement factors in sjn in the S/N 
due to filtering. 

Table IV. Accuracy of estimation of the parameters of the 2nd post- Newtonian signal for the three fiducial 
binaries. 



Binary 


S/N 




AM/M 




Am/m 


ASo/So 


ASs/Ss 


NS-NS 


15 


32 


0.080% 


42% 


63% 


60 X 10^% 


12 X 10*^% 


NS-BH 


32 


15 


0.26% 


40% 


59% 


11% 


19 X 10^% 


BH-BH 


77 


6 


0.92% 


150% 


230% 


240% 


890% 



We see that only the rms error in the chirp mass is small and also the accuracy of the determination of 
the spin-orbit parameter for NS-BH binary is satisfactory. The errors in reduced and total masses are 
large. 

One can derive simple general formulae for the accuracy of determination of the chirp mass, the reduced 
mass, and the total mass in terms of rms errors of the mass parameters ki. From the definition of the 
chirp mass one immediately obtains the following formula for the relative rms error in terms of the rms 
error in the mass parameter k, 

(55) 



AM/M = \rwoMpcAkM%'' . 



For the errors in the reduced and the total mass we obtain the following general formulae using the 
standard law of propagation of errors 



det 



Afi = 



Am 



ill / 



(56) 
(57) 



dk dki 
det = — — — 
m Ofj, 



(58) 



where 

dk dki 
dfj, dm 

and Ak, Aki arc rms error in mass parameters k and ki respectively. The formula above is the same when 
the 1st post- Newtonian, the 3/2 post-Newtonian, and the 2nd post-Newtonian corrections are included. 
We observe that errors in /x and m depend only on the masses and the rms errors in the parameters k 
and ki. The other mass parameters influence the errors in /i and m only through their correlations with 
the mass parameters k and ki and only through the functional form of the corrections as the rms error 
in the mass parameters are independent of their values. The errors in n and m are independent of the 
numerical values of the parameters k^^2 ^nd fc2- Since in general the rms error Afc is considerably smaller 
than Aki we get the following simplified expressions for the relative errors in the reduced and the total 
mass. 



A/i/ij, = ^riooMpcl^QAki, 
3 

Am/m = — riooMpcMoAfci, 



(59) 
(60) 



where a = 743/336 - 33/8 /it/m. We see that the error in the determination of the reduced mass and the 
total mass is determined by error in the first post-Newtonian mass parameter ki. Since the ratio fi/m is 
< 1/4 to a fairly good approximation we can take the value of a roughly equal to 1. 

If the spin effects could entirely be neglected and we would only have the reduced mass and the total 
mass as unknown in the mass parameters ki then we could achieve the accuracies in the parameters of 
the signal summarized in Table V. We considered three fiducial binary systems and 2nd post-Newtonian 
signal but with spin-orbit and spin-spin parameters removed. Thus the number of parameters estimated 
is 2 less than for the signal considered in Table IV. 



Table V Accuracy of estimation of the parameters for 3 fiducial binary systems and the 2nd post- 
Newtonian signal but with spin parameters removed. 
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Binary 


S/N 


AtcUis 




A/i//i 


Am/m 


NS-NS 


15 


0.47 


0.82 


0.29% 


0.43% 


NS-BH 


32 


0.32 


0.47 


0.19% 


0.28% 


BH-BH 


77 


0.18 


0.24 


0.27% 


0.37% 



We see that if spin parameters could be neglected we would have an excellent accuracy of estimation of 
the reduced and the total mass of the binary. 



2.4 The effects of other theories of gravity 

We shall consider two alternative theories. One is Jordan-Fiertz-Brans-Dicke (JFBD) theory (see for 
a detailed discussion) and the other is a multi-scalar field theory recently proposed in p^ . 

In the JFBD theory in addition to the tensor gravitational field there is also a scalar field. The theory can 
be characterized by a coupling constant that we denote by iv. General relativity is obtained when ui goes 
to infinity. The JFBD theory has two effects on gravitational emission. It admits dipole gravitational 
radiation and secondly there is a modification of the quadrupole emission due to the interaction of the 
scalar field with gravitating bodies. In the case of binary system the effects of the JFBD theory has 
been studied in great detail and a general formula for the change of orbital period was derived 
([^ eq. (14.22)). From that formula we get the following expression for the characteristic time r of the 
evolution of the binary due to radiation reaction in the case of circularized orbits and assuming that the 
contribution due to the dipole term is small 



(61) 



- 5 1 g^/' 1 
^ ^ 96 ^m2/3 K (7r/)8/3 ^ 

^ 192 K (7rm/)2/3^' 

where 

kB = (62) 

g = l-^(Ci+C2-CiC2), (63) 

- = g\l-^+'^l% (64) 

miC2 + m2Ci 

7 = 1 : , 65 

mi + 1712 

E = C1-C2. (66) 

Ci and C2 are "sensitivities" of the two bodies to changes of the scalar field. For a black hole the sensitivity 
C is always equal to 1. For a neutron star C depends on the equation of state. For neutron stars the 
sensitivity has been studied in ||37| for a number of equations of state and it was found for a wide range 
of such equations that it is proportional to the mass of the neutron star with proportionality constant 
varying from .17 to .31. Here we shall assume that Ci — 0.21miQ for a neutron star of m^Q solar mases. 
From the above formulae one sees that the dipole radiation will vanish if the binary system consists of 
two black holes or the neutron stars in the binary are the same. 

The Fourier transform of the signal in the stationary phase approximation including contributions due 
to JFKB theory is given by (we neglect any contributions due to eccentricity) 

h = Af-'^/^expi[2TTfta- <j)-TT/4: (67) 

5 

+ ^(a(/; fa)k' + ai(/; fa)h + a3/2(/; 70)^3/2 + a2(/; fa)ad{f; fa))kd)] 

for / > and by the complex conjugate of the above expression for / < where the function ad{f; fa) 
due to dipole radiation has the form 

(r 5 9 1 3 / _ 3 1 

ad(J,Ja) 192^0 (7r/)7/3 + 10 (^/„)io/3 7(^/,)7/3^' 
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and where 



1 Q^l"^ 

k' - 571^, (69) 

1 a^/^ 9 

Curent observational tests constrain lo to be greater than 600 and from timing of binary pulsar a lower 
limit on uj of 200 can be set. Thus it is sufficient to keep only the first terms in l/w. Then the two 
parameters above are approximately given by 

k' = -^{l-dkd), (71) 

kd = -^kB^' (72) 

where ^ 

dka^kB[^{C,+C2-C,C2) + ^-^]. (73) 

We thus see the the JFBD theory introduces a new parameter kd due to the dipole radiation and modifies 
the standard chirp mass parameter k by fraction dk^- 

We have investigated the potential accuracy of estimation of the parameter k^ assuming that the spin 
effects are negligible. We have taken neutron star/black hole binary with parameters given in Table I at 
the distance of 200Mpc. The result is summarized in Table VI. 

Table VI. The rms error for signal parameters in JFKB theory assuming spins are negligible for the 
binary of 1.4 solar mass neutron star and 10 solar mass black hole. 



S/N 


Aic [nis] 


A0C 




Am/m 




32 


0.47 


0.97 


0.57% 


0.73% 


2.3 X 10-^ 



The potential accuracy of determination of the dipole radiation parameter kd is high. Current obser- 
vational constraints indicate however that this parameter is small. We have the following numerical 
values. 

We conclude that the gravitational-wave measurement by planned long arm laser interferometers have 
the potential of testing JFBD to the accuracy comparable to tests in solar system and measurements 
from the binary pulsars [^6j . 

From the general class of tensor-multi-scalar theories studied recently |^ we shall consider a two- 
parameter subclass of tensor-bi-scalar theories denoted by T(/3',/3"). Theories in this subclass have 
two scalar fields and they tend smoothly to general theory of relativity when both parameters /?' and (3" 
tend to zero. The subclass is defined in such a way that the dipole radiation vanishes. From the general 
formulae one can calculate the characteristic time r. For circularized orbits the only modification is 
an effective change of the chirp mass parameter k given by the following formula 



k" = k-doF, (76) 
5 

5 



'^DF = —Ko{mi,Ci,m2,C2) + -iKq{mi,Ci,m2,C2)+ (77) 



Kdi(mi,Ci,m2,C2)) -f —Kd2{mi,Ci, 1112,02) 
4o 
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where coefficients Kq, '^q, ^^dii '«d2 ai'e due to contributions from quadrupole helicity zero, corrections to 
quadrupole lielicity two, and dipole radiation respectively. They are compUcated functions of the masses 
and sensitivities. We give the detailed formulae in Appendix C. In all the tensor- multi-scalar field theories 
whenever one of the component is a black hole corrections to the radiation reaction vanish. We have 
also found that for a simple model where sensitivities are proportional to masses of neutron stars and 
the proportionality constant is the same the correction doF does not depend on the parameter (3" . For 
a system of two identical neutron stars the correction cLdf takes a simple form 

doF = 0.21/3C^ (78) 

where C is the sensitivity of the neutron star to changes of the scalar field introduced above. Current 
observations constrain parameter (3 to be less than 1. For circularized orbits (the case considered above) 
the bi-scalar theory does not introduce a new mass parameter in the phase of the signal but only a shift 
in the "Newtonian" mass parameter k. We shall consider the possibility of estimating this shift in the 
next section. 



3 Search templates 
3.1 The Newtonian filter 

We have seen in the previous Section that the accuracy of estimation of the parameters is significantly 
degraded with increasing number of corrections even though a correction may be small. If we include the 
2nd post-Newtonian correction and filter for all unknown parameters then the accuracy of determination 
of the masses of the binary becomes undesirably low. Moreover we cannot entirely exclude unpredicted 
small effects in the gravitational-wave emission (e.g. corrections to general theory of gravity) that we 
at present cannot model. Thus there is a need for simple filters or search templates that will enable 
us to scan the data effectively and isolate stretches of data where the signal is most likely to be |jll|. 
The simplest such filter is just a Newtonian waveform which Fourier transform in stationary phase 
approximation is given by 

= ^2^z^^—^r"' expz[2^/t, - 0, - 7r/4 + k^^{^f)-^'\ (79) 

We shall call the Newtonian filter the filter which Fourier transform is given by the above formula and 
we shall denote it by N f . This filter has been investigated by the present authors ||, |3^, ^ and also by 
other researchers |2[ 0|. A different search template based on the post-Newtonian signal has 

recently been introduced in Rcf . |l4| . We discuss this alternative search template in the next subsection. 

In this section we examine the performance of the Newtonian filter. We demonstrate that such a template 
will perform well in detecting the signal from a binary and it also gives a reasonable idea of the nature of 
the binary. We shall investigate the performance of the Newtonian filter both analytically and numerically. 

Let us consider the correlation of the post-Newtonian signal with the Newtonian filter. Such an integral 
has the same form as the correlation integral given by Eq.^ in Section 2.1 except that all post- Newtonian 
mass parameters will be unmatched by the parameters of the filter. The correlation will be high if we 
can reduce the oscillations due to the cosine function as much as possible. Since the integrand of the 
correlation integral is fairly sharply peaked (HWHM ~ 26Hz) around its maximum at the frequency 
~ 47Hz we can achieve this by making the phase as small as possible around the peak frequency f^. 
The argument $ of the cosine in the integrand of the correlation of the post-Newtonian signal with the 
Newtonian filter including the effects due to eccentricity and dipole radiation takes the form 

$(/) = 27r/At + A(j)+ (80) 
5 

-^Hf-Ja)^k + a^if:, fa)ke + ai{f; fa)ki + a^i2{f\ fa)h/2 + a2{f] fa)k2 + aa{f-Ja)kd]- 

First we note that for all the mass parameter frequency functions ai{f]fa) the functions and their first 
derivatives vanish at the frequency fa- We shall therefore choose /„ — f^. Let us also transform the 
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phase parameter according to transformation given by Eq,53 with = f^. In the new parametrization 
the phase <i> takes the form 

$(/) = 2^(/-/:)At' + A</>"+ (81) 
^[a(/; /^)AA; - a,{f; Qk, + ai(/; Qk, - a3/2(/; Oh/2 + a2(/; /:)fc2 - ad(/; 

Let us examine the functional behaviour of $(/) around the frequency f'^. We find 

<i>(/)~2^(/-/^)At' + A</)"+ (82) 

5, , -2r 15''' , fcl ^3/2 , ^2 5 fed ^ 

'\J/Jo ) 



ggW/.o ^ L(^//)5/3 24 (tt/^) 19/9 tt/^ (^/^)2/3 (^/^)l/3 192 (7r/^)7/3 J 

+ o[(///:-i)^]. 

We see that in the above approximation we can make the phase $ vanish to the order (///^ — l)'^ when 
the following conditions hold 

^tmax = tpmax — t' ~ 0, (83) 
= ^Fma. - 4>" = 0, (84) 
^kmax — kpuiax ^ k = (S^) 

+ - hM-f'o) + hi^O^'^ ~ il (^ 

where subscript Fmax means the value of the parameter of the Newtonian filter that maximizes the 
correlation. Hence we can expect to match the Newtonian template to the post-Newtonian signal with 
the Newtonian mass parameter k shifted from the true value by a certain well-defined amount. The 
shift depends both on the parameters of the two-body system and the noise in the detector through the 
frequency /q. However the value of the shift in the k parameter is independent of the choice of the time 
parameter and phase in the Newtonian filter. 

In the following table we have given the numerical values of the shift in the parameter k calculated from 
Eq.^ for the 3 binary systems considered in the previous section. We have given three values of the 
shifts including one ((5fci), two ((5fc3/2), and finally three {5k2) post-Newtonian corrections. 

Table VII. Numerical values of the shifts in the mass parameter of the Newtonian filter calculated from 
the analytic formula (Eq^t 



Binary 


5ki 


5k3/2 


Sk2 


NS-NS 
NS-BH 
BH-BH 


0.03328 
0.01641 
0.004660 


0.01512 
0.005052 
0.001276 


0.01597 
0.006023 
0.001775 



We have also investigated the problem numerically and we have found the maxima to be located at the 
values of the shifts in the phase, the time, and the mass parameter k given in Tables VIIIA (1st post- 
Newtonian shift), VIIIB (3/2 post-Newtonian shift), VIIIC (2nd post-Newtonian shift) below. We have 
also given the factor I which is defined as 



In a previous work by these authors (|^, |^) we have claimed the factor I to be the drop in the signal-to- 
noise ratio as a result of using non-optimal (Newtonian) filter. However the signal-to-noise ratio falls as 
square of the factor (see Appendix D). We also give the range of integration over which we calculated 
the correlation. We have found that the we gain very little by extending the integration beyond that 
range. For the case of a neutron star binary increasing the range of integration up to 800Hz increases 
the signal-to-noise ratio by less than 1%. The reason for this is the effective narrowing of the band of the 
detector by the chirp signal discussed in the previous section. 

^We are grateful to T. Apostolatos for pointing this to us. 
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Table VIII. Numerical values of the factor I and shifts in the parameters of the Newtonian filter with 
respect to the true values for various post-Newtonian orders calculated numerically by maximizing the 
correlation function. 



A 



Binary 


h 


Ski 


St' 


S^" 


Range 


NS-NS 


0.68 


0.03721 


3.0 X 10-^ 


0.61 


30Hz - 200Hz 


NS-BH 


0.76 


0.01867 


1.7 X 10^3 


-0.53 


30Hz - lOOHz 


BH-BH 


0.85 


0.004931 


-4.1 X 10-3 


-0.20 


30Hz - lOOHz 



B 



Binary 


h/2 


Sk3/2 


St' 


S(f>" 


Range 


NS-NS 


0.90 


0.01564 


-3.5 X 10-3 


-0.40 


30Hz - 200Hz 


NS-BH 


0.87 


0.004905 


0.61 X 10-3 


0.068 


30Hz - lOOHz 


BH-BH 


0.87 


-0.001219 


0.39 X 10-3 


0.030 


30Hz - lOOHz 


C 


Binary 


h 


Sk2 


St' 


S^" 


Range 


NS-NS 


0.85 


0.01658 


-5.5 X 10-3 


-0.44 


30Hz - 200Hz 


NS-BH 


0.87 


0.006014 


-1.1 X 10-3 


-0.024 


30Hz - lOOHz 


BH-BH 


0.87 


0.001789 


-0.51 X 10-3 


-0.018 


30Hz - lOOHz 



We see that the agreement between the predicted values of the shifts in the parameters and the numerical 
values given above is very good. In particular the difference between the predicted values and the values 
of the shifts for the k parameter obtained numerically differ by less then 5%. 

The results of the detailed analysis carried out in p| show that when the amplitude and phase modu- 
lations due to the time dependence of the spin parameters are taken into account then in the worst case 
I — 0.63 for the correlation of the Newtonian filter with the 3/2pN signal. 

We have also performed the correlation using the signal in the time domain and evaluating the correlation 
using the fast Fourier transform. We kept the amplitude Newtonian. As we have remarked earlier the 
restricted post-Newtonian approximation are not equivalent in the frequency and the time domain. So 
the results are not the same. 
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Table IX. Numerical values of the I factor and the shifts obtained from the correlation of the Newtonian 
template with the signal in the time domain at various post-Newtonian orders. 



Binary 


k 


5ki 


h/2 


(5^3/2 


k 


Sk2 


Range 


NS-NS 


0.67 


0.04097 


0.97 


0.01576 


0.88 


0.01899 


30Hz - 200Hz 


NS-BH 


0.87 


0.01916 


1.00 


0.004889 


0.93 


0.01097 


30Hz - lOOHz 


BH-BH 


0.97 


0.005130 


1.00 


0.001896 


0.94 


0.005874 


30Hz - lOOHz 



We therefore conclude that the Newtonian filter will perform reasonably well in detecting the post- 
Newtonian signal. 

Using the Newtonian filter we would not like to loose any signals. We can achieve this by suitably lowering 
the detection threshold when filtering the data with the Newtonian filter. By this procedure we would 
isolate stretches of data where correlation has crossed the lowered threshold. The reduced data would 
contain all the signals that would be detected with the optimal filter but would also contain false alarms 
which number would bo increased comparing to number of false alarms with the optimal filter. This is 
the effect of lowering the threshold. The next step would be to analyse the reduced set of data with more 
accurate templates and the initial threshold to make the final detection. 

In Table X we have given examples of the performance of the above procedure. We assume the signal-to- 
noisc ratio threshold = 5 and we assume wc have 1 signal for the optimal signal-to-noisc ratio d. N is 
the expected number of detected signals with the optimal filter, Np is the number of false alarms, Nn is 
the number of detected signals with the Newtonian filter, Tjv is the lowered threshold, Nl is the number 
of signals with the lowered threshold and Np l is the number of false alarms with the lowered threshold 
(see Appendix D for definition of these qualities). 
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Table X. Comparision of number of true events and false alarms obtained with the optimal filter and the 
Newtonian filter. 



d 


FF 


N 


Nf 


Nn 


Tn 


Nl 


Nfl 


15 


.81 


27 


0.055 


20 


4.5 


28 


0.16 


15 


.36 


27 


0.055 


5.6 


3.225 


28 


2.1 


30 


.81 


225 


1.1 


165 


4.5 


230 


2.2 


30 


.25 


225 


1.1 


31 


2.875 


229 


32 



The theory of filtering with a suboptimal filter is outhned in Appendix D and the terms used in this 
Section are precisely defined. 

We have also calculated the covariance matrix for the parameters estimated with the Newtonian filter. 
Calculating the second derivatives of the correlation function at the maximum given by the numerical 
values of the parameters in Table VIII one gets the F matrix. The inverse gives the covariance matrix. 
The square roots of a diagonal components of the covariance matrix give lower bounds on the accuracy 
of determination of parameters with the Newtonian filter and they are approximate rms error for high 
signal-to-noise ratio as explained in Section 2. The results are summarized in Table XI for our three 
binary systems located at the distance of 200Mpc. The numbers are given for signals with the currently 
known post- Newtonian corrections but without the eccentricity and the dipolc terms. 

Table XI. Accuracy of determination of parameters of the Newtonian filter for the three fiducial binaries 
located at the distance of 200Mpc. 



Binary 


AtaN [ms] 




NS-NS 
NS-BH 
BH-BH 


2.9 
0.53 
0.22 


0.37 X 10-^ 
0.051 X 10-^ 
0.021 X 10^3 



One can easily calculate from Table II that the accuracy of determination of the mass parameter k with 
the Newtonian filter lies between the accuracy of determination of k for 1 and 3/2 post-Newtonian signal. 

In Appendix E we have derived a useful formula for the correlation function based on the approximation 
to the phase <& considered above. 

We shall next show that the Newtonian filter can also give a useful estimator characterizing the binary 
system. From the analytic investigation of the Newtonian filter given above it is clear that we can obtain 
an estimator of an effective mass parameter ks of the binary system given approximately by (cf.Eq.^) 

and the numerical investigation has shown that the Newtonian filter will determine the effective mass 
parameter which numerically value is accurately given by the above analytic formula. The ks parameter 
can be used to give an estimate of the chirp mass of the binary system. We define generalized chirp mass 
Mg as 

Mg = 1/4/' (88) 

We have calculated numerically the generalized chirp mass using the analytic formula (|8^) and we have 
found that it deviates from the true value by less than 4% for the range of masses from 1.4 to 10 solar 
mass. For the range of masses from 1.01 to 1.64 which is the expected range of neutron star masses given 
present observations of binary pulsars |Q the generalized chirp mass is always less than the true one by 
around 4% but with a very small range of .5% around the average value. 

Because of the inequality m > 2^^^A4 and the closeness of the generalized chirp mass to the true chirp 
mass the generalized chirp mass A4g gives a lower bound on the total mass of the system. Thus from its 
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estimate we can determine what binary system we observe. Also the R.H.S. of the above inequaHty gives 
a poor man's estimate of the total mass. For the range of masses of (1Mq,10Mq) it deviates by 50% 
from the true value of the total mass but for the range of (I.OIM0, 1.64Mq) acceptable for neutron star 
binaries it is only 5% smaller than the true mass. 

Another application of this estimate is that it can be used as an additional check on whether we are 
observing the real signal. If our estimate would deviate unusually from the predicted range of A4g 
corresponding to the range of individual masses of {1Mq,10Mq) we could veto the detection. 

An interesting application of the Newtonian filter would be to determine unexpected effects in the binary 
interaction that we would not be able to model and introduce into multiparameter numerical templates 
because we do not know their form. The idea is to use the estimates of the effective mass parameter 
kE- Particularly useful would be estimates of fc^ in the case of neutron star binaries. Since the range 
of the neutron star masses in a binary system is rather narrow the range of the allowable values for the 
generalized chirp mass will also be narrow. From the analysis in |^5| the range from the least lower 
bound and to the greatest upper bound is (I.OIMq, I.64M0) and the range from greatest lower bound 
to least upper bound is as narrow as (1.34M0, 1.43Mq). This implies the respective ranges in fc^ to be 
(0.57, 1.26) and (0.79, 0.71). From the population of estimates of the parameter we can determine 
its probability distribution and also the mean, variance or range of observed values of fc^;. One can then 
compare the observed distribution of kE and its characteristics with the ones obtained from observations 
of the neutron star binaries in our Galaxy or from the theoretical analysis and search for differences. As 
an example we consider Damour-Esposito-Farese bi-scalar tensor theory described at the end of Section 
2.4. The shift in the Newtonian mass parameter k due to effects of this theory is given by formula ([TS]). 
We have calculated this shift numerically and we have found that for the range of neutron star masses 
(I.OIM0 , 1.64Mo) and the parameter (3 = 1 (current observational bound) the shift is in the range of 
(0.018, 0.022). This shift is much larger than rms error in estimation of kE of 0.00037 (see Table XI). 
Consequently the effects of the bi-scalar theory could be determined to an accuracy depending on how 
well we would know the probability distribution of the neutron star masses and the number of available 
detections of gravitational waves from binaries. 



3.2 Post-Newtonian search templates 



In a recent work Jl4| different search templates than the Newtonian filter were recommended and ex- 
tensively analysed. The proposed templates are the post-Newtonian waveforms with all the spin effects 
and parameters removed. They have four parameters: amplitude, phase, reduced mass, total mass. We 
shall denote such search templates by IPNf, 3/2PNf, 2PNf where the number in front refers to the order 
of post-Newtonian effects included. In Ref.|14| the fitting factor FF {FF = P see Appendix D) of the 
3/2PNf search template was calculated and it was concluded that this template family works quite well 
even for signals with with both spin-modulational and the nonmodulated 3/2 post-Newtonian effects 
combined. In this Appendix we investigate the performance of the 2PNf search template for the case 
of the 2nd post- Newtonian signal in the approximation considered in Section 2. This means that we 
ignore all post-Newtonian effects in the amplitudes of both the signal and the template and we assume 
that the spin-orbit and the spin-spin parameters So and Ss in the signal are constant. In Table XII we 
give the factor / and the shift in the time parameter, phase, reduced mass and total mass for the three 
representative binary systems described in Section 2. We have also given the shifts in the reduced and 
the total mass parameters in percentages of their true values. 



Table XII. Performance of the 2nd post- Newtonian search template for the three fiducial binaries located 
at the distance of 200Mpc. 



Binary 


1 






5m 


dm 
ni 


(5t[ms] 


5(1) 


NS-NS 


0.98 


0.0028 


0.5% 


-0.017 


0.61% 


-9.5x10"^ 


0.00027 


NS-BH 


0.95 


0.52 


42% 


-4.8 


42% 


-3.0 


-0.28 


BH-BH 


0.98 


1.9 


38% 


-7.8 


39% 


2.3x10-3 


-0.00053 



We see that the 2PNf search template fits the signal better than the Newtonian search template Nf 
investigated in Section 3.1. There are two reasons for this. The 2PNf template has one more parameter 
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than Nf template and the phase of 2PNf template has all post-Newtonian frequency evolution terms 
whereas the phase of the Nf template has only Newtonian frequency evolution /~^/^. Also in the case 
of NS-NS binary which has small spin parameters the expectation values of the estimates of the reduced 
and the total masses are close to their true values. 

The advantage of the Newtonian search template might be its simplicity: it has the least possible number 
of parameters and hence the least computational time is needed to implement such a template in data 
analysis algorithms. Before the detailed data analysis schemes are developed for the real detectors it is 
useful to investigate theoretically a wide range of possible search templates. 

We have also calculated the covariance matrix for the 2PNf template. The results are summarized in 
Table XIII where we have given the rms errors in the time, reduced mass and the total mass parameters 
of this search template for the three binary systems. We have also given the errors in the reduced and 
the total mass in percentage of their true values. 

Table XIII. The rms errors in the estimators of the parameters of the 2nd post- Newtonian search 
template for the three fiducial binary systems located at the distance of 200Mpc. 



Binary 


Ata [ms] 


A^ipn[Mq] 


AmpN[MQ] 




Ampjv 
m 


NS-NS 


0.80 


0.0078 


1.1% 


0.011 


0.39% 


NS-BH 


0.40 


0.012 


1.0% 


0.0068 


0.06% 


BH-BH 


0.16 


0.0090 


0.2% 


0.0050 


0.03% 



We see that the rms errors of the parameters of the post-Newtonian search template are comparable to 
rms errors obtained with optimal filtering of the signal with spin parameters removed. 



4 Conclusions 

The analysis of the accuracy of estimation of parameters of the 2nd post-Newtonian signal (Section 2.3) 
has shown that main characteristics of this signal: chirp mass and the time parameter can be estimated 
to a very good accuracy: chirp mass to 0.1% - 1.0% amd time parameter to a quarter of a millisecond 
for typical binaries. A typical binary consists of compact objects of 1.4 to 10 solar masses and is located 
at the distance of 200Mpc from Earth and the amplitude of its gravitational wave signal is averaged 
over all directions and orientations. The signal-to-noise ratio of typical binaries varies from 15 to 77 for 
the planned advanced LIGO interferometers. However the accuracy of determination of post-Newtonian 
effects is considerably degraded due to large number parameters: 6 parameters in the phase of the 2nd 
post-Newtonian signal (Table II). Consequently the errors in determination of the reduced mass and the 
total mass are large and range from 50% to 200% for typical systems (Table IV). If spin effects could be 
neglected thereby reducing the number of parameters by 2 the rms errors of estimation of reduced and 
total masses would have a very impressive value of a fraction of a percent (Table V). 

Analysis of the accuracy of estimation of the effects of the dipole radiation in the Jordan-Fiertz-Brans- 
Dicke theory of gravity has shown that the planned laser interferometric gravitational wave detectors 
should have ability of testing alternative theories of gravity comparable to that of current observations 
in the solar system and our Galaxy. 

The numerical analysis of Section 2 supports the need for the search templates emphasized in Ref . |]ll| . 
The results of Section 3 show that the Newtonian filter (a search template with only one mass parameter) 
will perform reasonably well at least for the case of of constant spin parameters. Such a filter can be used 
to perform an on line scan of the data to search for the candidates for real signals. The measurement of 
the mass parameter of the Newtonian signal provides an accurate estimate of an effective mass parameter 
kE of the binary (see Eq.^). The value of this parameter gives the information about the binary 
analogous to the chirp mass in the analysis of the signal in the quadrapole approximation. Moreover this 
parameter contains information about the post-Newtonian effects and it can contain information about 
the effects that we cannot at present model for example about the effects due to unknown corrections to 
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general relativity in the strong field regime. Such information can be extracted if we built a probability 
distribution of from its estimators by the Newtonian filter. The post-Newtonian search templates 
analysed in perform better than Newtonian filters and considering increasing computational capability 
they can also be used in the on line analysis of the data. In the case of large spin parameters it would 
be useful to obtain relations of the two mass parameters in such templates to the true masses and spins 
similar to relation of the effective mass parameter of the Newtonian filter to the other parameters of the 
binary (see Eq. ^). For the case of the observed binary systems, binaries consisting of two neutron 
stars with small spin parameters the Newtonian filter will provide an accurate estimate of the chirp mass 
whereas the post- Newtonian search templates will provide accurate estimates of reduced and total masses. 
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Appendix A: The effects of eccentricity 



In this appendix we derive the first order correction due to eccentricity in the phase of the gravitational 
wave signal from a binary system. The derivation is due to N. Wex [ p8[ . 

Let a and e be respectively the semi-major axis and the eccentricity of the Keplerian orbit of a binary. 
From the quadrupole formula one obtains the following expressions for the secular changes of a and e 
averaged over an orbit 



/da\ 



g 1 + ZSgZ + 37 4 



\dt/ a3 (1 - e2)V2 



de\ 304 /3 e (1 + jle^) 
dt/ 15 a4 (l-e2)5/2 



(89) 



where (i = ^ m? fx. From these equations we get da/de which can be integrated with respect to e. The 
result is: 



a(e) = ao 



e(e) 

^(eo)' 



e(e) 



(^ _L 121 ^2 \ 870/2299 
12/19 \^ 304 / 



l-e2 



(90) 



where eo is an arbitrary initial eccentricity and ao = a(eo). From Kepler's third law tt/ = m^/^a'^/^, 
where / is gravitational wave frequency we get an analytic expression for / as a function of e. 



/ N /'i , 121 ^2\ 1305/2299 

/(e) = /o .(e) . e-/-li±iMO 



'7(e) ' 



(1 - e2)3/2 



(91) 



where /o — /(eo). For small eccentricities we find 



eo 



-19/18 



[l + 0{el) 



(92) 



Thus to first order in e the quantity I^, = Bq/q^^^ is a constant. We call le the asymptotic eccentricity 
invariant. The characteristic time for the evolution of the binary system is given by 



dfldt 



V da c\i) 



(93) 



From Kepler's third law we find 



96 ^m2/3 (7r/)8/3 1 



(1 - e'fl^ 

, 73 2 I 37 4 
^ 24^ ^ 96^ 



(94) 



For small eccentricities e we get 



96 ^m2/3 (7r/)8/3 



(95) 



Therefore using Eq.^we can express the characteristic time with first order correction due to eccentricity 
as 



1 



1 



96 ^m2/3 (7r/)8/3 



1 



157 
'2A 



=0 



-19/9' 



(96) 



The phase of the Fourier transform of the signal in the stationary phase approximation is given by 



V[J] = 2^ ft, - - ^/4 - 277 / T,(/')(l - ///') d/ = 



/ 



= 27r/ta - - 



1 



2/3 



577/ 



785 
1462 



^2/ J. N19/9 

eo(7r/o) 



+ 



(^j)5/3 (^J^)8/3 (^J^)5/3 

3477/ 43 



(v7/)34/9 (7r/,)43/9 (^/,)34/9 



(97) 
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and consequently the Fourier transform of our signal in the stationary phase approximation has the form 



Hf) = A.r'/' expi 2TTfta - - 7r/4 + ^{a{f; fa)k + a,{f- fa)k,) 
(and by the complex conjugate of the above expression for / < 0) where 



A = 



k = 



1 



1 ///2mi/3 



301/2 7,2/3 

1 



R 



1 



a{f-Ja) 



2/3 ' 

9 



^^2/3 
3 TTf 



40 (7r/)5/3 + 8 (^/,)8/3 5(7r/J5/3^ 

9 ^/ 



,„ , . 157 ^ 81 1 

aeUUa) = 77T 



24 1^1462 (7r/)34/9 + 43 (7r/a)43/9 34 (7r/a)34/9 



for / > 0. (98) 



(99) 
(100) 
(101) 
(102) 



We have investigated the accuracy of measurements of parameters of the above signal with first order 
eccentricity contribution. We have considered neutron star/neutron star binary. The results 
rized in Table XIV. 



i are summa- 



Table XIV The rms errors of the parameters of the signal with first order contribution due to eccentricity 
for a binary of two neutron stars of 1.4 solar mass each at the distance of 200Mpc. 



S/N 


Ate [ms] 




A/i//i 


Am/m 


Afce[MQ^/^(100Hz)i9/9] 


15 


0.56 


1.2 


0.50% 


0.74% 


3.6 X 10-^ 



However for the currently observed binaries the eccentricity invariant is extremely small. For Hulse- 
Taylor pulsar Jg — 1.8 x 10~^'^[Mq^^^100}1z^^/^]. We have the following numerical values. 

k, = 1.3 X lo-"(j^^-iij-^)(iji)V3|„-V»io(,H.'»/»] (103) 

^.2.8.10W(T^^)(4^f)=", (104) 

where r2ooMpc is distance in 200Mpc. Thus for eccentricity effects to be measured one would need 
extremely short period binaries of high eccentricity. Such binaries could perhaps occur in the center of a 
galaxy or be created as a result of some supernova explosions. 



Appendix B: Covariance matrices at various post-Newtonian orders 



In this Appendix we give the numerical values of the covariance matrices at various post-Newtonian 
orders for the reference binary. The reference binary has the chirp mass of 1 solar mass and is located 
at the distance of lOOMpc. We only give reduced covariance matrices i.e. covariance matrices for the 
phase parameters. As indicated in Section 2 the estimator of the amplitude parameter is uncorrelated 
with phase parameters. The integration range in the Fisher matrix integrals was taken to be from lOHz to 
infinity and the spectral density of advance LIGO detectors was assumed (Eq.p5|). The frequency fa was 
chosen such that the rms error in the time parameter is minimum. The minimum frequency is denoted 
by frn and its numerical value is given for each covariance matrix. The subscripts N, IPN, 3/2PN, 2PN, 
2PNe refer to signal including quadrupole radiation, 1st post-Newtonian correction, 3/2 post-Newtonian 
correction, 2nd post-Newtonian correction, and 1st order effect due to eccentricity respectively. The order 
of parameters in the matrices is the following: tc, (t>c, k, fci, k^/2i ^2, kg. 

.C = 70Hz 



/ 1.96x10-^ 8.2x10-*^ 1.21 X IQ-i" \ 
Cat = 8.2 X 10-6 0.00537 2.08 x 10"^ (105) 
\ 1.21 X 10-1° 2.08 X 10-^ 6.65 x IQ-" / 
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PIPN _ 



lOOHz 



CiPN 



Jm = 160Hz 



C: 



3/2PJV 



/ 2.24 X 10-^ 9.93 x 10"^ 

9.93 X 0.00754 

8.87x10-1° -6.28x10-^ 

\ -6.61 X 10-*^ 0.000147 



3.87 X 10-^° 
-6.28 X 10"^ 
1.4 X 10-^ 
-1.97 X lO-'^ 



-6.61 X 10-8 \ 

0.000147 
-1.97 X 10-^ 

0.0000289 / 



3.41 X 10"** 

0.00002 
-9.65 X 10"^ 
3.36 X lO"'^ 
0.0000212 



0.00002 
0.0191 
-9.34 X 10"^ 
-0.000403 
-0.007 



-9.65 X 10"' 
-9.34 X 10" 

2.23 X 10"* 
-8.71 X 10-' 

-0.0000622 



3.36 X 10"'' 
-0.000403 
-8.71 X lO"'' 
0.00349 
0.0253 



(106) 



0.0000212 
-0.007 

-0.0000622 
0.0253 
0.185 



(107) 



p2PN 



lOOHz 



C2PN 



\ 



5.62 X 10"** 
0.0000286 

2.39 X lO-'' 
-9.48 X 10"' 
-0.000199 
-0.00102 



0.0000286 

0.0189 
0.000016 
-0.0162 
-0.26 
-1.11 



2.39 X 10"" 
0.000016 
2.3 X 10"^ 
-0.000162 
-0.00219 
-0.00796 



-9.48 X 10" 

-0.0162 
-0.000162 
0.116 
1.59 
5.86 



6 



-0.000199 

-0.26 
-0.00219 

1.59 

22. 

81.5 



-0.00102 

-1.11 
-0.00796 

5.86 

81.5 

305. 



(108) 



f2PNe ^ i20Hz 
C2PNe — 
I 



6.36 X 10"* 


0.0000369 


-3.56 X 10"* 


4.9 X lO"'^ 


-0.000135 


-0.000977 


-4.86 X 10"" 


0.0000369 


0.0279 


0.000035 


-0.0315 


-0.474 


-1.89 


-5.21 X 10"" 


-3.56 X 10"* 


0.000035 


2.36 X 10"** 


-0.00121 


-0.0143 


-0.0464 


1.33 X 10"'' 


4.9 X 10"^ 


-0.0315 


-0.00121 


0.63 


7.57 


24.7 


-6.52 X 10"^ 


-0.000135 


-0.474 


-0.0143 


7.57 


91.5 


301. 


-7.59 X 10"" 


-0.000977 


-1.89 


-0.0464 


24.7 


301. 


999. 


-0.000024 


-4.86 X 10"" 


-5.21 X 10"" 


1.33 X 10"" 


-6.52 X 10"^ 


-7.59 X 10"^ 


-0.000024 


8.28 X 10"" 



(109) 



Appendix C : Coefficients in the Damour-Esposito-Farese biscalar T(/3',/3") theory 

The coefficients Kq, k^, k^i, Krf2 in the shift of the Newtonian mass parameter k due to the biscalar 
T(/3',/3") theory (Eq.^ in Section 3.4) are given by the following formulae 

= \0B{Ci+Cl), (110) 
= p'B{Cfx2+Clxi), (111) 

Kdl = ^f3'B{C^,Xi-Clx2){xi-X2), (112) 

Kd2 = {abi2i - ab22i)xi + {ab2i2 - abii2)x2, (113) 

where 

XI = ^, (114) 
m 

X2 = ^, (115) 
m 

and constant A and B have the values 

A = 2.1569176, 5 = 1.0261529. (116) 
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Ci and C2 are sensitivities of the two bodies to changes of the scalar field. The functions ab are given by 



abi 



ab 



'212 



ab 



221 



(3'{-C2 - BCl + {A- 3B)Ci - (A - B)2C2C^ 



(117) 



(2^2 - TAB + 5B^)CiCf) + (3'^B^{-3Cf + 2CfCi +Cf + ^CzCf + AC^Ct) + ^P" BC 



P'i-Ci - BCl + iB)Cl -{A- B)2CiC 



{2A^ - TAB + 5B^)CfCi) + fi'^B'^[-3Cl + 2CfCi + + ^CiC| + ACfCl) + ^l3"BCf 



(118) 



/3'{-C2 - \b{CI + C|) + {A- 3B)Cl -{A- B)C2{Cf + C|)+, 



2 



{2A' - TAB + 5B')Ci{Ct + Ci)) + (3' B\-3C'i + Ci{Ci + Ci) 

2r'i\ I -'^ o" D/-'2 



2 ^ ^C|C? + AC?C4) + ]^p"BCl 
1 



a6n2 = (i'{^C2 - ^B{Cl + Cl) + (A - 3B)Cf - (A - B)Ci(C? + C|)+, 

i(2A2 - 7^5 + hB^)Cl{Cl + CD) + /3'^b2(-3Cj* + Cl{Cl + Cl) + 

Ct + \clCl + AClCt) + \p"BCl. 
For a detailed exposition of the theory the reader should consult Ref. |^ . 

Appendix D: Detection of the known signal v^rith a non-optimal filter 



(119) 



(120) 



Suppose that we would like to know whether or not in a given data set x there is present a signal h. We 
assume that the noise n in the data is additive. There are two alternatives: 



NO SIGNAL 
SIGNAL 



X — n 

X = h + n 



(121) 



A standard method to determine which of the two alternatives holds is to perform the Neyman-Pearson 
test[^8| . This test consists in comparing the likelihood ratio A, the ratio of probability density distributions 
of the data x when the signal is present and when the signal is absent, with a threshold. The threshold 
is determined by the false alarm probability that we can tolerate (the false alarm probability is the 
probability of saying that the signal is present when there is no signal). The test is optimal in the 
sense that it maximizes the probability of detection of the signal. In the case of Gaussian noise and 
deterministic signal h the logarithm of A is given by 



InA = {x\h) - ^{h\h). 



(122) 



Thus in this case the optimal test consists of correlating the data with the expected signal and it is 
equivalent to comparing the correlation G := {x\h) with a threshold. The probability distributions po 
and pi of G when respectively the signal is absent and present are given by 

Po{G;d) = ^=L=exp ^ 



PiiG;d) 



%/27rd2 
1 



: exp 



{G-d^ 



d^ 



(123) 
(124) 



V2^ 

where d is the optimal signal-to-noise ratio d^ — {h\h) and we assumed that the noise is a zero mean 
Gaussian process. 

Let T be a given threshold. This means that we say that the signal is present in a given data set if G > T. 
The probabilities Pp and Pd of false alarm and detection respectively are given by 



PF{T,d) 
PDiT,d) 



PoiG;d) dG, 
Pi{G;d) dG 



(125) 
(126) 
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In the Gaussian case they can be expressed in terms of the error functions. 

PpidT,d)^^erM^), 
P^(ciT,rf) = i(l + er/(^^), 



(127) 
(128) 



where erf and erfc are error and complementary error functions respectively [ [47| . and we have introduced 
for convenience the quantity dx '■= VT that we call the threshold signal-to-noise ratio. In practice we 
adopt a certain value of the false alarm probability that we can accept and from formula (127) we calculate 
the detection threshold T. 

Let _F be a linear filter and let n be the additive noise in data x then 



{x\F)^is\F) + in\F). 



The signal-to-noise (S/N) ratio is defined by 



(S/Nf := 



E,[is\F)'] is\Ff 



E,[{n\Fr] {F\F) 



(129) 



(130) 



where Ei means expectation value when the signal is present. By Schwartz inequality we immediately 
see that (S/N) is maximal and equal to d when the linear filter is matched to the signal i.e. F = h. This 
is another interpretation of the matched filter - it maximizes the signal-to-noise ratio over all linear filters 
|p8| . However when the noise is not Gaussian the matched filter is not the optimal filter; it does not 
maximize probability of detection of the signal. We see that in the case of Gaussian noise the problem of 
detecting a known signal by optimal filter is determined by one parameter - the optimal signal-to-noise 
ratio d. 

Suppose that because of certain restrictions of practical nature we cannot afford to use the optimal filter 
h and we use a suboptimal one - which is not perfectly matched to the signal. Thus (/i|ft,jv) < {h\h). 
We denote y^7i|7ijv) by do and we assume that (/iAr|ft.Ar) — {h\h) — d^. Our suboptimal correlation 
function is given by Gn — (a;|/ijv) and its probability distributions pNo and pNi when respectively the 
signal is absent and present are given by 



PmiGN^d) 
PNi{GN\d, do) 



1 



^/2^ 
1 



exp 



exp 



{Gn 



j2\2 



d^ 



(131) 

(132) 



We see that the suboptimal detection problem is determined by two parameters - d and do, square roots 
of the expectation values of the optimal and suboptimal correlations when the signal is present. The false 
alarm and detection probabilities as in the optimal case can be expressed in terms of the error functions. 



PpidT,d)^^erfci^), 
P^[dT,d,do) = \{l + erf(^^). 



(133) 
(134) 



We see that the probability of false alarm for the suboptimal case is the same as in the optimal case 
however the probability of detection in the suboptimal case is always less then the probability of detection 
in the optimal case since do < d and the error function erf(x) is an increasing function of the argument 
X. The signal-to-noise ratio in the case of suboptimal linear filter h^ is given by 



(S/Nf 



{h\hNf 

(/ijvl^w) 



(135) 



Let us denote the ratio do/d by I. The ratio I measures the drop in the expectation value of the correlation 
function as a result of non-optimal filtering. We see that due to suboptimal filtering the signal-to-noise 



29 



ratio decreases by square of the factor /. We denote P by FF and following Rcf. call it the fitting 
factor. 

In our considerations we need to calculate the number of events that will be detected by linear filtering. 
We shall make a number of simplifying assumptions. We shall assume a Euclidean universe where in the 
sphere of radius Tq we have one source and that at the distance the optimal signal-to-noise ratio is d. 
Moreover we shall assume that the magnitudes of the signal h and the suboptimal filter hjsr are inversely 
proportional to the distance r from the source. Then the square roots dr and dor of the expectation 
values of the optimal and suboptimal correlations at the distance r are given by 

dr = —d, (136) 
r 

T 

dor = —do- (137) 
r 

We assume that the sources are uniformly distributed in space. Then the expected number of detected 
real events N and Nj\j in the optimal and the suboptimal case respectively is given by 



47r r^PD{dT,dr)dr 



3 o 



/•oo 

3 / x'^PD{dT,d/x)dx (138) 
Jo 

47r r^PNoidr, d, do/r) dr 



3 o 



3 / x'^PND{dT,d,do/x)dr. (139) 



NidT,d) 



NNidT,d,do) = 





The assumptions that led to the above formulae mean that we neglect general relativistic, cosmological 
and evolutionary effects. Because of the noise even if there is no signal there is always a non zero 
probability that the correlation function crosses the threshold. Thus there will be a certain number Np 
of false events. For a given optimal signal-to-noise ratio and a threshold dx this number is the same for 
both the optimal and suboptimal filter and it is given by 

47r f°° r^PF(dT,d/r)dr 

- " 4^ 3 = (140) 

3 o 



/>oo 

3 / x^PpidTjd/x) dx. 

J o 



We observe that in the Gaussian case the integrals in the formulae (|13§| ), (139), and (141) are convergent 
even though we integrate over the all infinite Euclidean volume. 

Appendix E: An approximate formula for the correlation function. 

In this Appendix we shall derive an approximate formula for the correlation integral. Let us consider 
the expression for the correlation function given by (^6|) . The integrand of the correlation integral is the 
product of the integrand of the signal-to-noise integral Ind{f) considered in Section 2 and oscillating 
factor. We know that the Ind{f) is a fairly sharply peaked around a certain frequency consequently 
to obtain a reasonable approximation we expand the phase around the frequency /q. Keeping only the 
terms to the second order we get 

$(/) ^ 27r(/ - QM' - A0"-t-, (141) 

where 

C = - 2TrfX (142) 
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and 

We shall call Ue an effective mass parameter. Aks is the difference in the effective mass parameter of the 
signal and the filter. Thus in the above approximation the post-Newtonian signal can be parametrized 
by one effective mass parameter ks- In other words the dimension of the parameter space of the filters 
is effectively reduced. This last interpretation has been emphasized in where 1st post-Newtonian 
corrections to the phase were considered. The mass parameter estimated by Newtonian filter considered 
in Section 3 is just the effective mass parameter. We stress that the parameter depends not only on 
the parameters of the two-body system but also on the characteristic frequency of the noise in the 
detector. 

The next step is to obtain a manageable approximation to the function Ind{f). We approximate it by 
a Gaussian function with the mean equal to the frequency and the standard deviation equal to the 
HWHM (Tq of the function Ind{f). We extend the range of integration from —oo to -|-oo. We introduce 
a normalization factor such that the integral of the approximate integrand is equal to the optimal signal- 
to-noise ratio d. It is then useful to introduce a reduced correlation integral H' — H/d^ where d is the 
S/N ratio. Thus our approximate formula for the reduced correlation integral takes the form 

K = -— / exph(/ - /:)/(2af )] cos[27r(/ - fJAt' + A0" -I- ^llIf>AkTB] (144) 

The above integral can be done analytically. It is convenient to introduce the following new variables 
and new parameters 

V = (145) 



= A0", (146) 



2^At'^2a',\ (147) 
' (148) 



96/f(^/^)^/3 
then our integral takes a simple form 

1 /■+°° 

H'^{d, t,k)^— / exp[-?/2] cos[-d" + ry + Ky^]dy. (149) 

J-oo 

We see that in the new variables introduced above the reduced correlation integral is independent of the 
characteristics of the integrand Ind{f ) i.e. and <j'^. The analytic formula for the function H'^({}" , r, k) 
is given by 

Kir, r, n) ^ (YT^ '"P[-4(lT^] cos[l/2(arctan. - ^(^^) " (150) 

By appropriate transformations given in Section 2 we can obtain approximate formulae to the corre- 
lation integral for an arbitrary choice of the time and the phase parameters. Let us first consider the 
transformation given by Eq.(53). In the coordinates introduced above it takes the form 

r^d- l=Tp (151) 

Then the approximate formula for the correlation function is given by 

r, .) = ^^^1^ -P["i(T^] cos[l/2(arctan. - ^^^) + i=rp - ^] (152) 

We see that in these new coordinates for k = the correlation function oscillates with the maxima at the 
discrete values of r coordinate given by 

Tmax = i^y (153) 
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where 

P = /o/^o (154) 
and n is an integer. In the original coordinates Eg. ( 153 ) takes the form Ai — 1//q. Thus the correlation 



integral oscillates with the period determined by the characteristic frequency of the noise of the detector 
f 

The expressions for the correlation function for different choice of the time and the phase parameters can 
be obtained by the following transformations. These are transformations given by Eqs.(^9|) and expressed 
in our dimensionless coordinates. 

^ = l?' + -«;y(l-(55/3), (155) 

5 

t = t' + ^kV(1-(5«/3) (156) 

where 

5 = f'olfa- (157) 

From the approximate formula for the correlation function obtained above we see that the correlation 
is given by the product of an oscillating cosine function and an envelope. In the cosine function there 
are oscillations with the period of 1//^. The envelope function is exponentially damped if we move away 
from the maximum at the center except for the direction given by r = along which the damping is 
least. The equation of the ridge r = in the primed coordinates is given by 

^ k'p'{1-5^'^), (158) 



4V2 

and in the original coordinates it takes the form 

Consequently we conclude that the general appearance of the correlation function in coordinates At' and 



Afc is a series of peaks aligned along a straight line given by Eq.(159) above and occurring with the 
period l//o in the time coordinate. Numerical investigation shows that the correlation integral exhibits 
these properties and that our analytic formula reproduces qualitatively its behaviour. The approximate 
formula obtained above may be a useful tool for developing algorithms to recognize the chirp signal in a 
noisy data set. 
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